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We performed angle resolved photoemission (ARPES) studies on Cu(lll) and Au(lll) surface 
states with circularly polarized light. Existence of local orbital angular momentum (0AM) is con- 
firmed as has been predicted to be broadly present in a system with an inversion symmetry breaking 
(ISB). The single band of Cu(lll) surface states is found to have chiral 0AM in spite of very small 
spin-orbit coupling (SOC) in Cu, which is consistent with theoretical prediction. As for Au(lll), 
we observe split bands for which 0AM for the inner and outer bands are parallel, unlike the Bi2Se3 
case. We also performed first principles calculation and the results are found to be consistent with 
the experimental results. Moreover, majority of 0AM is found to be from d-orbitals and a small 
contribution has p-orbital origin which is anti-aligned to the spins. We derive an effective Hamil- 
tonian that incorporates the role of 0AM and used it to extract the 0AM and spin structures of 
surface states with various SOC strength. We discuss the evolution of angular momentum structures 
from pure 0AM case to a strongly spin-orbit entangled state. We predict that the transition occurs 
through reversal of 0AM direction at a /c-point in the inner band if the system has a proper SOC 
strength. 
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I. INTRODUCTION 

When solids possess both inversion and time reversal 
symmetries, Kramer's theorem dictates that each elec- 
tronic state is doubly spin degenerate.^ When the inver- 
sion symmetry is broken at surfaces and interfaces, one 
can expect that the spin degeneracy is lifted except at 
some special /c-space points. Lifting of spin degeneracy 
and resulting band splitting in two dimensional electron 
gas such as metallic surface states are typically explained 
in terms of the Rahsba effect.^ In the Rashba effect, elec- 
tron spin interacts with an effective magnetic field stem- 
ming from electron motion in a surface electric field and 
results in Zeeman splitting. Consequently, the surface 
state band splits and attains a chiral spin structure at 
a constant energy contour. Surface band splitting and 
concomitant chiral spin structure have been experimen- 
tally observed on surfaces of relative high atomic num- 
ber metals^^, interfaces of hetero-structures^ and surface 
states of topological insulators.^ 

In spite of its success in explaining the energy split- 
ting and chiral spin structure, the original Rashba model 
could not provide the proper energy scale, giving about 
10^ times smaller value than the measured one in the 
case of Au(lll) surface bands. There have been sev- 
eral theoretical studies to resolve the issue but they did 
not address all the aspects of the Rashba effect j2rii It 
was only recently found that local orbital angular mo- 
mentum (GAM) plays the key role in the Rashba- type 
band splitting by inducing asymmetric charge distribu- 
tion in Bloch states when the atomic spin-orbit coupling 
(SOC) is much larger than the crystal field energy The 
asymmetric charge distribution interacts with the surface 



electric field and provides the energy scale. It results in 
chiral GAM structure and chiral spin structure naturally 
follows from the strong SGC. A surprise came when chi- 
ral GAM was found to exist even if there is no SGC.^« In 
this case, contrary to the strong SGC case, GAM vectors 
for the degenerate state are found to be parallel to each 
other while spins are anti-parallel. 

As discussed above, the spin and GAM configurations 
and band splitting are quite different for the weak and 
strong SGC cases. Then, it would be interesting to in- 
vestigate the transition from weak SGC to strong SGC 
cases. It is especially meaningful because the most stud- 
ied Rashba split surface states on Au(lll) possibly be- 
long to an intermediate SGC case. To address this issue, 
we performed angle resolved photoemission (ARPES) ex- 
periments on Cu(lll) and Au(lll) surface states with 
circularly polarized light as well as first-principles calcu- 
lation. We also develop an effective Hamiltonian to study 
the problem in an analytical way. We confirm that GAM 
is indeed present in the Cu(lll) and Au(lll) surface 
bands with most contribution from d-orbitals. Analysis 
based on the effective Hamiltonian shows that transition 
from parallel GAM in the weak SGC case to anti-parallel 
GAM in the strong SGC case occurs through reversal of 
GAM at a /c-point in the inner surface band. 



II. METHODS 

ARPES measurements were performed at the beam 
line 9A of HiSGR equipped with VG-SCIENTA R4000 
analyzer. Data were taken with right and left circularly 
polarized (RCP and LCP) 10 eV photons. The total en- 
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ergy resolution was set to be 10 meV at 10 eV, and the 
angular resolution was 0.1°. We performed experiment at 
10 K under a base pressure better than 7.5 x 10~^^ Torr. 
To obtain clean and well-ordered surfaces, we cleaned the 
surfaces by using Ar sputtering and annealed the sam- 
ples by e-beam heating. In order to check high qual- 
ity of the sample surface, we performed low energy elec- 
tron diffraction (LEED) and confirmed a long range or- 
der. For the density-functional theory (DFT) calcula- 
tions within the local-density approximation (LDA), we 
used the OpenMX^^ based on the linear-combination-of- 
pseudo-atomic-orbitals (LCPAO) method. Spin-orbit 
interaction was included via the norm-conserving, fully 
relativistic j-dependent pseudopotential scheme in the 
non-collinear DFT formalism. To calculate the 0AM 
and spin angular momentum (SAM) for a specific k- 
point, we used the LCAO coefficients of local atoms. 
Due to the non-orthogonality of pseudo-atomic-orbitals, 
LCAO coefficients are not strictly normalized to unity. 
To compensate for this, we re- normalized the coefficients 
by assuming the orthogonality, and obtained SAM and 
0AM values that are strictly bounded above by 0.5 and 
1, respectively. 



III. RESULTS AND DISCUSSION 

In figure [TJ we plot ARPES results from Cu(lll) sur- 
face states. As expected, Cu(lll) surface states have spin 
degeneracy and show a single parabolic dispersive band 
due to the small atomic SOC. The binding energy at 
the point is E(r) = 418 meV, and the Fermi momentum 
kp = 0.21174"-*^ (see figure. 1(a) - 1(d)). These values are 
consistent with published values. In figure 1(e), we plot 
the circular dichroism (CD) at the Fermi energy defined 
as CD = RCP - LCP. The CD data presented in the 
color scale changes gradually from red to blue as value 
increases continuously from — to +. Plotted in figure 
1(f) is the cut image along ky = shown as the dashed 
line in figure 1(e). This figure shows that CD is nega- 
tive (positive) for positive (negative) ky for all energies. 
The band split near the point is an artifact resulted from 
broadening in the LCP data due to aging. It is clear from 
the raw data in figures 1(b) and 1(d) that there is only a 
single band. Finally, plotted in figure 1(g) is the normal- 
ized CD (defined as NCD = (RCP - LCP) / (RCP + 
LCP)) at constant binding energies as a function of the 
azimuthal angle defined in the inset. The curves have a 
sine function form, which suggests that the 0AM forms 
a chiral structure. ^^^^^ In addition, the estimated magni- 
tude of 0AM from CD decreases as the binding energy 
increases. The magnitude is found to be approximately 
proportional to the magnitude of the electron momentum 
value. 

The behavior of CD in ARPES reveals that chiral 
0AM indeed exists in the Cu(lll) surface band in spite 
that SOC is very small in Cu as predicted. It is also con- 
sistent with the prediction that the magnitude of 0AM is 
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FIG. 1: ARPES results from Cu(lll) surface states, (a) 
Fermi surface and (b) the cut along the kx = line (dashed 
line in panel (a)) taken with RCP light, (c) and (d) The same 
for LCP light, (e) and (f) RCP - LCP data, (g) Normalized 
CD defined as NCD = (RCP - LCP) / (RCP + LCP) as a 
function of the azimuthal angle defined in the inset. 



linear in the electron momentum value k. The formation 
of chiral 0AM is a way of lowering the system energy by 
making the charge distribution asymmetric in the pres- 
ence of surface electric field. As will be discussed later, 
0AM vectors for the degenerate band are parallel to each 
other while spins are anti-parallel. Note that if CD were 
due to spins, we would have not observed CD because 
there is no spin polarization for a state. The fact that 
we can observe CD from Cu(lll) shows that CD is from 
OAM^^ 

We now turn our attention to the Au(lll) surface 
states, possibly the most studied surface states in re- 
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FIG. 2: ARPES results from Au(lll) surface states, (a) 
Fermi surface and (b) the cut along the kx = line (dashed 
line in panel (a)) taken with RCP light, (c) and (d) The same 
for LCP light, (e) and (f) RCP - LCP data, (g) NCD as a 
function of the azimuthal angle. Data from the two bands are 
summed in the estimation to have averaged NCD. 



gards to the Rashba effect. The first direct experimental 
evidence for Rashba split bands was from the Au(lll) 
surface.^ This system has a similar band shape to that of 
Cu(lll) surface states, but has band splitting of about 
110 meV. In figure 2(a), we plot Fermi surface map of 
Au(lll) surface states taken with RCP light. Clear dou- 
ble Fermi surfaces due to Rashba-type splitting are seen 
in the figure. Figure 2(b) shows the cut along the kx = 
(dashed line in figure 2(a)). Once again, clear split bands 
are seen. The observed band bottom is at E(r) = 479 
meV and the Fermi vectors are kp = 0.165A~^ for the in- 
ner band and kp = 0.196A~^ for the outer band. These 



values are quite consistent with reported values. Data 
taken with LCP light in figures 2(c) and 2(d) show sim- 
ilar features compared to the RCP data except there is 
some difference in the intensity profile. 

Figures 2(e) and 2(f) plot circular dichroism RCP — 
LCP for the Fermi surfaces and bands along the kx = 
line (dotted line in 2(e)). Except that there is clear band 
splitting, the overall CD profile looks similar to that of 
Cu(lll) : it is negative (positive) for ky>0{<0) region. 
Averaged NCD of the inner and outer bands as a function 
of the azimuthal angle plotted in figure 2(g) also shows 
a similar behavior to that of Cu(lll) that it can be fit 
with a sine function and the magnitude is approximately 
linear in k. However, we also note that NCD value for 
Au is about three times larger than that of Cu, which 
suggests that 0AM in Au case is generally larger than 
that for Cu. A peculiar aspect to note is that, unlike 
the Bi2Se3 case,^^ inner and outer bands have the same 
CD sign even though their spin directions are opposite. 
Same sign of CD suggests that 0AM of the inner and 
outer bands are pointing in the same direction as in the 
Cu(lll) case in spite of a stronger SOC and thus larger 
splitting. On the other hand, a careful look reveals that 
CD for the inner band is slightly stronger by about 35% 
compared to that from the outer band. 

In order to investigate the spin and 0AM structures in 
more detail, we performed first principles density func- 
tional theory (DFT) calculations for Au(lll) surface 
states within the local density approximation. In the 
DFT results plotted in figure 3, top panels (3(a) and 
3(b)) show directions and sizes of spin and 0AM by the 
arrows. The left panel 3(a) is for the inner band while the 
right one 3(b) is for the outer band. The spin and 0AM 
directions obtained from DFT results are as expected. 
The directions for 0AM (blue) and spin (red) are op- 
posite to each other in the inner band while the outer 
band has parallel 0AM and spin. The spin and 0AM 
directions are consistent with previous report^ and ex- 
perimental result in figure 2, respectively. That is, while 
spins are anti-parallel between the two bands, 0AM are 
pointing in the same direction. 

We turn our attention to the magnitudes of 0AM and 
spin. We first look at the spin magnitude. It is seen from 
the figures that the spin magnitudes from the two bands 
are quite similar (only the directions are opposite). More- 
over, it has very little momentum dependence. In fact, 
the magnitude actually slightly decreases as the momen- 
tum increases. As for 0AM, we find that 0AM magni- 
tude increases as we move away from the F point, making 
0AM magnitude approximately linear in k as indicated 
in the experimental results. However, the 0AM magni- 
tudes for the inner and outer bands are also very similar, 
which appears to contradict the experimental result. 

To better understand the seemingly contradicting re- 
sults from experiment and theory, we look into atomic 
orbital dependent contributions to 0AM. For Au(lll) 
surface states, 5d- and 6p-orbitals contribute to OAM^^. 
In figures 3(c) and 3(d), we plot contributions from 5d- 
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FIG. 3: DFT results on the angular momentum structures 
of Au(lll) surface states. DFT results for (a) inner and (b) 
outer bands. 0AM (spin) of a state is represented by blue 
(red) arrow. Shown in the lower panels are contributions from 
p-state (red arrow) and d-state (blue dashed arrow) for (c) 
inner and (d) outer bands. 



and 6p-orbitals (blue and red, respectively). One can see 
that the d-orbital contribution dominates and determines 
the 0AM direction. On the other hand, p-orbital contri- 
bution shows quite different behavior from that of d. The 
direction is opposite for inner and outer bands, making 
it always anti-parallel to the spin direction. Even though 
the total 0AM is similar for the inner and outer bands, 
they are built differently and result in different circular 
dichroism. As a side note, the composition of the 0AM 
shown in figures 3(c) and 3(d) suggests that the atomic 
SOC parameter for 6p is larger than that of hd because a 
large value of atomic SOC parameter tends to anti-align 
the spin and 0AM. 

A natural question is how the spin and 0AM structures 
evolve as a function of SOC strength, from a small SOC 
(e.g., Cu) to strong SOC (e.g. Bi2Se3). To elucidate the 
issue, we wish to develop an effective Hamiltonian and 
analyze the evolution. It has already been shown that 
free electron based model cannot explain various aspects 
of split bands and that tight binding state is needed.^ A 
through derivation with nearest neighbor hoppings con- 
sidered on tight binding states can be found elsewhere^^. 
Instead, we wish to develop a simpler effective model for 
the surface states. For simplicity, we limit our discussion 
to the p-orbital case, but it can be extended to other 
orbit als. 

For electrons in the surface states, there are 4 terms 
that are significant in the Hamiltonian. They are the 
kinetic energy Rk^ atomic spin orbit coupling HsoCi 
crystal field HcF-, and the electrostatic energy due to 
interaction of asymmetric charge distribution with sur- 



face electric field Hes-^^ For the kinetic energy Hk^ we 
simply add a /c^ term at the end to account for the free 
electron-like parabolic band. Hsoc is (^L • S where a 
is the atomic SOC parameter. Meanwhile Hcf is A 
for px and p^, and for py (note that we take y-ax.is 
as surface normal, not the usual z-axis). The last term 
Hes comes from interaction of the surface electric field 
and electric dipole moment of a state with asymmetric 
charge distribution.^^ The asymmetric charge distribu- 
tion is a combined effect of electron momentum and lo- 
cal 0AM, and should be proportional to the momentum 
and 0AM. Therefore, the Hamiltonian is given 
Hes = -P- Es -a^^L x k) - = -aj,{k x E^) • L 
where is a constant that is related to how efficiently 
asymmetric charge is created by k and L. It is similar to 
the well known Rashba Hamiltonian Hr = a^{k x Eg) • a 
with the spin operator replaced by the 0AM operator, 
but can account for the split energy. 

We consider a state along the x-direction k = kxX with- 
out loss of generality, in which case Hes — ~(^xkx 
The total Hamiltonian is estimated in the basis of \pxt)^ 
\Pxi), \Pyt)i \Pyl)i \Pzt)i and |p^;). The result is given 
by 

H = Ck^i^ 

/AO ia^kxEs-f 

A ia^kxEg + 




H,C. 

V 

The 6x6 matrix is diagonalized to obtain the eigen- 
states and energies. The two lowest energy states cor- 
respond to the surface states we have discussed above, 
and their spin and 0AM can be easily estimated. 

The spin and 0AM of the two states are obtained as 
a function of atomic SOC parameter a to investigate the 
evolution of spin and 0AM. The results are plotted in 
figure 4. Red (blue) color denotes spin (0AM) while 
dotted (solid) line is for the inner (outer) band, a in- 
creases as we go from left to right. For the small a case 
in figure 4(a), the band is degenerate. The 0AM and 
spin structures are consistent with the results discussed 
above. 0AM forms a chiral structure with the magnitude 
approximately linear in k while spins are anti-parallel to 
each other but co-linear with 0AM. This case represents 
the Cu(lll) surface states. We also look at the other ex- 
treme of the large a case in figure 4(c). In this case, the 
degeneracy is lifted with a large band splitting. 0AM 
and spin are always anti-parallel to each other due to the 
large a and the magnitude of 0AM is independent of 
the electron momentum /c, which is consistent with the 
results from Bi2Se3.~ 

We finally examine the intermediate case in figure 4(b) 
which should provide us information on how parallel 
0AM for the small a case evolves to anti-parallel 0AM 
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other hand, for a large k value, Hes is dominant and the 
system lowers the energy by inducing parallel 0AM as in 
the Cu(lll) case. In the case of Au, such reversal of the 
0AM in the inner band was not observed, which means 
a for Au may not be large enough. Observation of such 
0AM reversal may be possible for surface states of Pb. 



Increasing SOC parameter a 

FIG. 4: Spin and GAM structures calculated by using the 
effective Hamiltonian for (a) small, (b) intermediate and (c) 
large atomic SGG parameter a. All the parameters are fixed 
except a. We chose 0.195 as the Fermi momentum kp oi 
the outer band, similar to the case of Au(lll). The red and 
blue arrows represent spin and GAM, respectively. Dotted 
(solid) arrows are for the inner (outer) band. 

for the large a case. For the outer band, 0AM and spin 
are anti-parallel and the 0AM magnitude has a linear 
k dependence. This is similar to the small a case. The 
inner band, however, shows a quite different behavior. 
While 0AM and spin are parallel to each other near the 
Fermi energy, they become anti-parallel near the F point 
(small k value). When k increases from 0, 0AM gradu- 
ally decreases, reverses the direction and increases again. 
Therefore, the 0AM structure evolves from parallel con- 
figuration in the small a case to anti-parallel configura- 
tion in the large a case by reversing the 0AM direction 
in the inner band at a certain k point instead of reversing 
the direction at all k points. Such behavior stems from 
the fact that Hsoc is larger than Hes for a small k and 
0AM prefers to stay anti-parallel to the spin. On the 



IV. SUMMARY 

ARPES studies on Cu(lll) and Au(lll) surface states 
with circularly polarized light as well as first principles 
calculation have been performed to investigate the spin 
and 0AM structures. Experimental and theoretical re- 
sults show that 0AM indeed exists even for Cu(lll) and 
Au(lll) surface states as predicted earlier. Cu(lll) has 
a degenerate single band with chiral and parallel 0AM 
while Au(lll) has split bands with 0AM structure not 
too much different from that of Cu(lll). DFT calcula- 
tion shows that majority of 0AM comes from d-states 
while a small contribution is from p-orbitals. We also de- 
veloped an effective Hamiltonian with the role of 0AM 
incorporated to investigate the evolution of the spin and 
0AM structures as a function of the atomic SOC param- 
eter. We find that there should be 0AM reversal at a 
specific momentum in the inner band when the system 
has a proper spin-orbit coupling strength. 

This work is supported by NRF (contract No. 
20090080739) and the KICOS under Grant No. 
K20602000008. This work was performed with the ap- 
proval of the Proposal Assessing Committee of HSRC 
(Proposal No.lO-A-57). 



* Electronic address: changyoung@yonsei . a c . kr | 

^ H.A. Kramers, Proc. Acad. Sci. Amsterdam 33, 959 (1930). 

^ Y. A. Bychkov and E. I. Rashba, JETP Lett. 39, 78 (1984). 

^ S. LaShell, B. A. McDougall, and E. Jensen, Phys. Rev. 
Lett. 77, 3419 (1996). 

^ M. Hoesch, M. Muntwiler, V. N. Petrov, M. Hengs- 
berger, L. Patthey, M. Shi, M. Falub, T. Greber, and 
J.Gsterwalder Phys. Rev. B 69, 241401 (2004). 

^ Christian R. Ast and Hartmut Hochst, Phys. Rev. Lett. 
87, 177602 (2001); Yu. M. Koroteev, G. Bihlmayer, J. E. 
Gayone, E. V. Ghulkov, S. Bliigel, P. M. Echenique, and 
Ph. Hofmann, Phys. Rev. Lett. 93, 046403 (2004); T. Hira- 
hara, K. Miyamoto, 1. Matsuda, T. Kadono, A. Kimura, T. 
Nagao, G. Bihlmayer, E. V. Ghulkov, S. Qiao, K. Shimada, 
H. Namatame, M. Taniguchi, and S. Hasegawa, Phys. Rev. 
B 76, 153305 (2007). 

^ D. Pacile, G. R. Ast, M. Papagno, G. Da Silva, L. Mores- 
chini, M. Falub, Ari P. Seitsonen, and M. Grioni, Phys. 
Rev. B 73, 245429 (2006); Christian R. Ast, Jurgen Henk, 
Arthur Ernst, Luca Moreschini, Mihaela C. Falub, Daniela 
Pacile, Patrick Bruno, Klaus Kern, and Marco Grioni, 
Phys. Rev. Lett. 98, 186807 (2007); Christian R. Ast, 
Daniela Pacile, Luca Moreschini, Mihaela C. Falub, Marco 



Papagno, Klaus Kern, and Marco Grioni, J. Henk, A. 
Ernst, S. Ostanin, P. Bruno , Phys. Rev. B 77, 081407(R) 
(2008); Fabian Meier, Vladimir Petrov, Sebastian Guer- 
rero, Christopher Mudry, Luc Patthey, Jiirg Osterwalder, 
and J. H. Dil, Phys. Rev. B 79, 241408(R) (2009). 

^ Junsaku Nitta, Tatsushi Akazaki, Hideaki Takayanagi, and 
Takatomo Enoki, Phys. Rev. Lett. 78, 1335 (1997) 

^ M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 
(2010). 

^ L. Petersen and P. Hedegard, Surf. Sci. 459, 49 (2000). 
Miki Nagano, Ayaka Kodama, T Shishidou and T Oguchi, 
J. Phys. Condens. Matter 21, 064239 (2009). 
Koichiro Yaji, Yoshiyuki Ohtsubo, Shinichiro Hatta, Hi- 
roshi Okuyama, Koji Miyamoto, Taichi Okuda, Akio 
Kimura, Hirofumi Namatame, Masaki Taniguchi and Tet- 
suya Aruga, Nature Commun. 1, 1 (2010). 
^2 G. Bihlmayer, S. Bliigel, and E. V. Ghulkov, Phys. Rev. B 
75, 195414 (2007). 

E. Frantzeskakis, S. Pons, and M. Grioni, Phys. Rev. B 
82, 085440 (2010). 

Tamio Oguchi and Tatsuya Shishidou, J. Phys. Condens. 
Matter 21, 092001 (2009). 

Seung Ryong Park, Choong H. Kim, Jaejun Yu, Jung Hoon 



6 



Han, and Changyoung Kim, Phys. Rev. Lett. 107, 156803 
(2011). 

Jin-Hong Park, Choong H. Kim, Jun Won Rhim, and Jung 

Hoon Han, arXiv:1112.1821 

http://www.openmx-square.org. 

T.Ozaki, Phys. Rev. B 67, 155108 (2003). 

G. Theurich and N. A. Hill, Phys. Rev. B 64, 073106 

(2001). 

M. Mulazzi, G. Rossi, J. Braun, J. Minar, H. Ebert, G. 
Panaccione, I. Vobornik, and J. Fujii, Phys. Rev. B 79, 
165421 (2009). 



Seung Ryong Park, Jinhee Han, Chul Kim, Yoon Young 
Koh, Changyoung Kim, Hyungjun Lee, Hyoung Joon 
Choi, Jung Hoon Han, Kyung Dong Lee, Nam Jung 
Hur, Masashi Arita, Kenya Shimada, Hirofumi Namatame, 
Masaki Taniguchi, arXiv: 1103.0805 

Y. H. Wang, D. Hsieh, D. Pilon, L. Fu, D. R. Gardner, Y. 
S. Lee, and N. Gedik. Phys. Rev. Lett. 107, 207602 (2011) 
J. Henk, M. Hoesch, J. Osterwalder, A. Ernst and P. 
Bruno., J. Phys. Condens. Matter 16 7581 (2004) 



